Anomalous net-baryon-rapidity spectra at RHIC 
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Net-baryon rapidity distributions in central Au+Au collisions at ^snn = 200 GeV are 
shown to consist of three components. The nonequilibrium contributions are accounted for 
in a Relativistic Diffusion Model. Near midrapidity, a third fraction containing Z eq ~ 22 
protons reaches local statistical equilibrium in a discontinuous transition. It may be 
associated with a deconfmement of the participant partons and thus, serve as a signature 
for Quark-Gluon Matter formation. 
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Rapidity distributions of participant (net) baryons are very sensitive to the dynamical 
and statistical properties of nucleus-nucleus collisions at high energies. Recent results 
for net-proton rapidity spectra in central Au+Au collisions at the highest RHIC energy 
of yj SjviV — 200 GeV show an unexpectedly large rapidity density at midrapidity. The 
BRAHMS collaboration finds pQ dN/dy = 7.1 ± 0.7 (stat.) ± 1.1 (sys.) at y = 0. 

The A, A feed-down corrections reduce this yield by 17.5 per cent [Tj when performed 
in accordance with the PHENIX A— results j2| at 130 GeV, but the amount of stopping 
remains significant, although a factor of about 4 smaller as compared to Pb+Pb at the 
highest SPS energy. (A corresponding STAR result jS] for y=0 at 130 GeV does not yet 
include the feed-down correction). Many of the available numerical microscopic models 
encounter difficulties to predict the net-proton yield in the central midrapidity valley of 
the distribution, together with the broad peaks at the detected positions. 

In this Letter I propose to interpret the data in a nonequilibrium-statistical Relativistic 
Diffusion Model. The net baryon rapidity distribution at RHIC energies emerges from a 
superposition of the beam-like nonequilibrium components that are broadened in rapidity 
space through diffusion due to soft (hadronic, low p±) collisions and particle creations, 
and a statistical equilibrium (thermal) component at midrapidity that arises from hard 
(partonic, high p±) processes. 

At RHIC energies, the underlying distribution functions turn out to be fairly well 
separated in rapidity space. Since the transverse degrees of freedom are in (or very close 
to) thermal equilibrium, they are expected to decouple from the longitudinal ones. The 
time evolution of the distribution functions is then governed by a Fokker-Planck equation 
(FPE) in rapidity space [U El El HI IH1 Ej- In the more general case of nonextensive 
(non- additive) statistics [HI] that accounts for long-range interactions and violations of 
Boltzmann's Stosszahlansatz j^l ^] as well as for non-Markovian memory (strong 
coupling) effects JHE], the resulting FPE for the rapidity y in the center-of-mass frame 
is 



Since the norm of the rapidity distribution has to be conserved, p = 1 is required, and the 



0_ 

Of 



[R(y,t)} 




J(y)[R(y,t)r\+D(t) — [R(y,t)r. 



(1) 



2 



nonextensivity parameter that governs the shape of the power-law equilibrium distribution 
becomes q = 2 - v ^U|. In statistical equilibrium, transverse mass spectra and transverse 
momentum fluctuations in relativistic systems at SPS-energies v /J/v r /v=17.3 GeV require 
values of q very slightly above one, typically q = 1.038 for produced pions in Pb+Pb [10]. 
For q — > 1, the equilibrium distribution converges to the exponential Boltzmann form, 
whereas for larger values of q (with q < 1.5) significantly broader equilibrium distributions 
are obtained, and the time evolution towards them becomes superdiffusive [TTH I13j. 

To study rapidity distributions in multiparticle systems at RHIC energies in a 
nonequilibrium-statistical framework [U IHj , I start with q = v — 1 corresponding 
to the standard FPE. For a linear drift function 

J(y) = (y eq -y)/r y (2) 

with the rapidity relaxation time r y , this is the so-called Uhlenbeck-Ornstein process, 
applied to the relativistic invariant rapidity. The equilibrium value is y eq = in the 
center-of-mass for symmetric systems, whereas y eq is calculated from the given masses 
and momenta for asymmetric systems. Using 5— function initial conditions at the beam 
rapidities ±yb (±5.36 at p=100 GeV/c per nucleon), the equation has analytical Gaussian 
solutions. The mean values shift in time towards the equilibrium value according to 

< ViAt) >= W 1 - exp(-2t/r y )] ± y b exp (-t/r y ). (3) 

For a constant diffusion coefficient D y , the variances of both distributions have the well- 
known simple form 

<Tl 2 (t) = D v T v [l-exp(-2t/T v )], (4) 

whereas for a time dependent diffusion coefficient D y (t) that accounts for collective (mul- 
tiparticle) and memory effects the analytical expression for the variances becomes more 
involved [IT] . At short times t/r y « 1, a statistical description is of limited validity due 
to the small number of interactions. A kinematical cutoff prevents the diffusion into the 
unphysical region \y\ > y^. For larger values of t/r y , the system comes closer to statistical 
equilibrium such that the FPE is valid. 

Since the equation is linear, a superposition of the distribution functions emerging 
from Ri t2 (y, t = 0) = 5(y =F y^) yields the exact solution, with the normalization given by 
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the total number of net baryons and the value of t/r y at the interaction time t=Ti nt (the 
final time in the integration of (1)) determined by the peak positions 4J. This approach 
has also been applied sucessfully to produced particles at RHIC energies although 
there the initial conditions are less straightforward. 

The microscopic physics is contained in the diffusion coefficient. Macroscopically, the 
transport coefficients are related to each other through the dissipation-fluctuation theorem 
(Einstein relation) with the equilibrium temperature T 

D y = a-T^f( Ty ,T). (5) 

In ji] I have obtained the analytical result for D y as function of T and r y from the 
condition that the stationary solution of (1) is equated with a Gaussian approximation to 
the thermal equilibrium distribution in y-space (which is not exactly Gaussian, but very 
close to it) as 



D y (r y , T) — [c(V5, T)m 2 T ■ (1 + 2- + 2(-) 2 )} *ezp(™) (6) 
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with c(y/s,T) given in |5- in closed form. This allows to maintain the linearity of the 
model and hence, to solve the FPE analytically, although small corrections are to be 
expected. They cause minor deviations in the calculated rapidity distributions that are 
within the size of the error bars of the experimental data at SPS energies. 

In the linear model, net baryon rapidity spectra at low SIS-energies are well repro- 
duced, whereas at AGS-, SPS- and RHIC energies I find discrepancies to the data that 
rise strongly with y/s. At SPS energies, this has recently been confirmed in a numerical 
calculation (HI Ej based on a nonlinear drift 

J(y) = —a ■ m±sinh(y) = —a ■ p\\ (7) 



with the transverse mass m± = \Jm 2 + p 2 ±) and the longitudinal momentum Together 
with the dissipation-fluctuation theorem (5), this yields exactly the Boltzmann distribu- 
tion as the stationary solution of (1) for v — q — 1. The corresponding numerical solution 
with 5— function initial conditions at the beam rapidities is, however, only approximately 
correct since the superposition principle is not strictly valid for a nonlinear drift. Still, the 
numerical result shows almost the same large discrepancy between data and theoretical 
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rapidity distribution as the linear model. In a q=l framework, the net proton distribution 



in Pb+Pb at the highest SPS energy requires a rapidity width coefficient D y r y that is 
enhanced beyond the theoretical value (5) by a factor of g{\fs) — 2.6 due to memory and 
collective effects HEIDI], Fig.l. 

Alternatively, a transition to nonextensive statistics EH EI] maintaining the weak- 
coupling diffusion coefficient from (5) requires a value of q that is significantly larger than 
one. In an approximate numerical solution of (I) with the nonlinear drift (7), q=1.25 has 
been determined for the net-proton rapidity distribution in Pb+Pb collisions at the SPS 
jHllZl- The only free parameter is q, whereas in the linear q=l case the enhancement of 
D y beyond (5) is the only parameter. 

This value of q in the nonlinear model is considerably larger than the result 
g(y / sjvjv)=1.12 extrapolated from Wilk et al. jHj at the SPS-energy ^/snn = 17.3 
GeV. Here, the relativistic diffusion approach is applied to produced particles in proton- 
antiproton collisions in the energy range a/s = 53 GeV - 1800 GeV, and used to predict 
LHC-results. The nonlinearity q > 1 appears to be an essential feature of the pp data. 
The larger value of q in heavy systems as compared to pp at the same NN-center-of-mass 
energy emphasizes the increasing superdiffusive effect of multiparticle collisions both be- 
tween participants, and between participants and produced particles. It is, however, 
conceivable that both a violation of (5) due to memory effects, and q > 1 have to be 
considered in a complete description. 

Investigating the Au+Au system at RHIC energies in the linear q=l model, it is now 
easy to show that in central collisions (10 per cent of the cross-section) the nonequilibrium 
net-proton rapidity spectrum calculated with the standard dissipation-fluctuation theorem 
underpredicts the widths of the nonequilibrium fractions of the experimental distribution 
[TH significantly. In this comparison, the temperature T=170 MeV is taken from a 
thermal fit of charged antiparticle-to-particle ratios in the Au+Au system at 200 GeV 
per nucleon [TH] ITH] , and the theoretical value of the rapidity width coefficient calculated 



from the analytical expression (6) is J D y r y = 7.6 • I0~ 2 . As was shown in |6j for SPS 

results, the discrepancy persists in case of the nonlinear drift (7) that yields the exact 
Boltzmann-Gibbs equilibrium solution for q=l 
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An enhancement of the weak-coupling rapidity width coefficient by a factor of g(y/s) ~ 3.7 
due to collective and memory effects in the system corresponding to a violation of (5) 
yields a good reproduction of the nonequilibrium contributions with T int /T y =0.26, but 
the midrapidity valley that is present in the data is completely absent in the extensive 
nonequilibrium q=l case, solid curves in Fig.l (bottom). This remains true in the nonex- 
tensive case (1 < g < 1.5), with an approximate distribution function jTUl EH El IE] that 
is given by a linear superposition of Tsallis-like power-law solutions of (1) 

R h2 (y,t) = [1 - (1 - q)^cosh(y- < Vl , 2 (t) (9) 

The dashed curves in Fig.l show the result for q=1.4, T=170 MeV and a mean transverse 
mass < m±_ >=1.2 GeV. This solution is far from the nonextensive equilibrium distri- 
bution which would be reached for < yi :2 (t oo) >= y eq , and it is significantly below 
the midrapidity data. The result is even worse for larger values of m±. In contrast, the 
Pb+Pb data at SPS energies |201 are well described both in the linear model [3] (Fig.l, 
top) and in the nonlinear case (cf. 6J for results with a time-dependent temperature and 
an integration over transverse mass). 

It turns out, however, that the RHIC data can be interpreted rather precisely in the 
linear q=l framework with the conjecture that a fraction of Z eq ~ 22 net protons near 
midrapidity reaches local statistical equilibrium in the longitudinal degrees of freedom. 
The variance of the equilibrium distribution R eq {y) at midrapidity is broadened as com- 
pared to the Boltzmann result (dashed curve in Fig. 2) due to collective (multiparticle) 
effects by the same factor that enhances the theoretical weak-coupling diffusion coef- 
ficient derived from (5). This may correspond to a longitudinal expansion velocity of 
the locally equilibrated subsystem as accounted for in hydrodynamical descriptions. In 
the nonextensive model, the corresponding equilibrium distribution is broadened (blue- 
shifted) according to q ~ 1.4. 

Microscopically, the baryon transport over 4-5 units of rapidity to the equilibrated 
midrapidity region is not only due to hard processes acting on single valence (di) quarks 
that are described by perturbative QCD, since this yields insufficient stopping. Instead, 
additional processes such as the nonperturbative gluon junction mechanism [TJ] are nec- 
essary to produce the observed central valley. This may lead to substantial stopping 
even at LHC energies where the separation of nonequilibrium and equilibrium net baryon 
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fractions in rapidity space is expected to be even better than at RHIC. In the late ther- 
malization stage [TH], nonperturbative approaches to QCD thermodynamics are expected 
to be important. 

Recent work indicates that one may account for the observed stopping in heavy- 
ion collisions at SPS and RHIC energies with string-model parameters determined from 
hadron-hadron collisions ^H]- If this was confirmed, the corresponding rapidity distribu- 
tions would not be considered to be anomalous from a microscopic point of view. However, 
this view does not offer a distinction between nonequilibrium and equilibrium contribu- 
tions to the net baryon rapidity spectra, which both exist at RHIC energies, and are 
anomalously broadened. 

Macroscopically, the complete solution of (1) in the q=l case is a linear superposition 
of nonequilibrium and equilibrium distributions (Fig.2, bottom) 

R(y, t = T int ) = Rl(y, Tint) + R 2 {y, T int ) + Req{v)- (10) 

It yields a good representation of the preliminary BRAHMS data 1 . (In the q > 1 case, 
the corresponding solution is questionable because the superposition principle is violated). 
Based on (10), the transition from net-proton rapidity spectra with a central plateau in 
Pb+Pb at the lower SPS energies [TT], via a double-humped distribution at the maximum 
SPS energy [201 EJ 13 HI] to the central valley at RHIC [T] is well understood. It has not 
yet been possible to identify a locally equilibrated subsystem of net baryons at midrapid- 
ity below RHIC-energies, although it cannot be excluded that it exists. At SPS energies, 
the data |20] are well described by the nonequilibrium distributions, and it is much more 
difficult (and probably impossible) to identify a locally equilibrated component because 
the relevant rapidity region is comparatively small, and an equilibrated contribution can- 
not be separated from the nonequilibrium components in rapidity space. In pp-collisions 
at y/s =53-900 GeV, no convincing signatures of a phase transition were found [21] . 

Most remarkably, Fig.2 suggests that in central Au+Au collisions at x /Tnn= 200 GeV 
there is no continuous transition from the nonequilibrium to the equilibrium contribution 
in net-proton rapidity spectra as function of time. This may well be due to a sudden 
enhancement in the number of degrees of freedom as encountered in the deconfinement 
of participant partons, which enforces a very rapid local equilibration in a fraction of the 
system. The central valley in net-proton rapidity spectra at RHIC energies could thus be 
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used as an indicator for partonic processes that lead to a baryon transfer over more than 
4 units of rapidity, and for quark-gluon plasma formation. 

To conclude, I have interpreted recent results for central Au+Au collisions at RHIC 
energies in a Relativistic Diffusion Model (RDM) for multiparticle interactions based on 
the interplay of nonequilibrium and equilibrium ("thermal") solutions. In the linear ver- 
sion of the model, analytical results for the rapidity distribution of net protons in central 
collisions have been obtained. The anomalous enhancement of the diffusion in rapid- 
ity space as compared to the expectation from the weak-coupling dissipation-fluctuation 
theorem due to high-energy multiparticle effects has been discussed using extensive and 
nonextensive statistics. 

A significant fraction of about 14 per cent of the net protons reaches local statistical 
equilibrium in a fast and discontinuous transition which is likely to indicate parton de- 
confinement. The precise amount of protons in equilibrium is related to the experimental 
value of the rapidity density close to y=0 and hence, possible changes in the final data 
will affect the percentage. It has not yet been possible to isolate a corresponding fraction 
of longitudinally equilibrated net protons in the Pb+Pb system at SPS energies. Since no 
signatures of a transition to the quark-gluon plasma have been observed in pp-collisions, 
Quark-Matter formation is clearly a genuine many-body effect occuring only in heavy 
systems at sufficiently high energy density. Consequently, a detailed investigation of the 
flat midrapidity valley found at RHIC, and of its energy dependence is very promising. 
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Figure captions 



FIG. 1: Nonequilibrium contributions to the net-proton rapidity spectra of Au+Au at 
i/sjvat = 200 GeV in the Relativistic Diffusion Model (RDM) with an equilibrium 
temperature of T=170 MeV (bottom). The solid curve is obtained in the linear 
model (q=l) with an enhanced width E] due to multiparticle effects according 
to the preliminary BRAHMS data [I], squares. The dashed curve corresponds to 
q=1.4 and < m±_ >=1.2 GeV in the nonlinear model. At SPS energies (top), 
NA49 data [20J for central events (5 per cent) including A feed-down corrections are 
compared with the linear model [H lllj. 

FIG. 2: Net-proton rapidity spectra for central collisions of Au+Au at ^s^n — 200 
GeV consist of two nonequilibrium components (solid peaks, top) plus an equilib- 
rium contribution at T=170 MeV, dashed curve. It is broadened due to collective 
(multiparticle) effects, shaded area, and after hadronization, it contains Z eq ~ 22 
protons. Superposition creates the flat valley near midrapidity (bottom) in agree- 
ment with the preliminary BRAHMS data points pQ; diamonds include A feed-down 
corrections at y=0 (17.5 per cent) and y=2.9 (20 per cent), respectively. Arrows 
indicate the beam rapidities ±y b . 
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Fig.l. 
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Fig. 2. 



